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SUMMARY 



The analysis of various types of shell under combined 
"bending and torsion is discussed. The calculation and the 
use of the shear center are touched upon as incidental 
problems. Twelve fully worked numerical examples are give 
in an appendix. 



INTRODUCTION 



The literature on shell analysis is quite scattered 
and some of it is not easily available. A definite heed 
therefore exists for setting down in a reasonably compre- 
hensive, but concise, manner the principles and tho meth- 
ods used in the various phases of shell analysis. The 
present paper deals with the distribution of the stresses, 
chiefly the shearing stresses, over the cross sections of 
cantilever shells of constant cross section subjected to 
combined bending and torsion. 

The subject matter is far from being new; as the ap- 
pended bibliography indicates, the main principles were 
well established in 1931. Continued discussion in tho 
tochnical literature indicatos, however, that the knowl- 
edge of the subject is not 30 widely disseminated as it 
needs to bo. It is hopod that tho manner of presentation 
choson for this paper, in particular tho collection of 
fully worked numerical examples, will help materially to 
achieve the ultimate goal, namely, to provide all prac- 
ticing engineers with a working knowlodgo of the subject. 
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THE ANALYSIS OF SHELLS 



Basic Assumptions and Theories 



Typical cases of the problem, to be treated, are shown 
in figure 1(a). As indicated "by the figure, the shell 
will be assumed to have a constant crops section. It will 
also "be assumed, in general, that the material effective 
in "bending is disposed symmetrically about the horizontal 
axis. The horizontal and the vertical ares will then be 
principal axes, and the load will be assumed to act verti- 



For purposes of stress analysis, the structures are 
idealized in the usual manner. A certain effective width 
of- skin is added to each. actual longitudinal or flange to 
obtain the cross-sectional area effective in bending; this 
effective area is assumed fro be concentrated at the con- > 
troid. The skin itself is assumed to carry only shear. 
Idealized sections aro ropresented as in figure l(b). 

Stresses caus-ed by bending are obtained in first ap- 
proximation by applying the engineering theory of bending 
to the idealized crpss sections. This theory is based on 
the assumptions that plane cross sections remain plane and 
that Hooke's law applies. The theory leads to the formu- 
las, fcr the normal stresses due to bending, 



and, for the shear stresses due to bending, 



The derivation of these formulas can be found in any text- 
book on strength of mate-rials. 

It should be noted that, in the computation of the 
static moments Q, as well as of the moments of— dnertia I, 
consideration is given only to the material assumed to be 
effective in bending. 



cally . 




(2 a) 



In shell structures, it is often convenient to use, 
not the shear stress T, but the shear force per inch 
length of shoet, which will be designated by 
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q. - T t 

and which will "be called the "shear-force intensity" o 
"briefly, "shear intensity." Formula (2a) then becomes 



■q = 



for open sections,, where "b = t. 





-(a) 



■o 




tt>) . 

Figure 1. 
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In the case of a closed cross section, an equivalent for- 
mula is obtained "by considering the equilibrium of hori- 
zontal forces, on a piece of the crps„e_ section as shown in 
figure 2 



a - cl = (2h) 
^■m r r 



where Q, is now the static moment a"bout the neutral axis 
of the (effective) areas of the longitudinals lying "be- 
tween the skin panels r 'and m "where the shear intensi- 
ties are measured. 




Figure 2. 



Shear stresses caused by torsion in a tube (fig-r 3) 
are obtained by the well-known formula 

applicable to thin-wall tubes. The angle of twist per 
unit length of tube is given by 

9 = GJ 

where the torsion constant J is defined by 
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In practical cases, the thickness is constant- over large 
parts of the circumference; the calculation of the line 
integral therefore reduces to the addition of a few. frac- 
tions of the type s/t. Substitute (3) into the formula 
for angle of twist, and there is obtained 

0 2&A ft 



In the most general case, Gr and q. may be variable 
along the circumference. Variation of G- may "be due to 
the use of different materials or to the formation of 
diagonal-tension fields. Variation of q. may be caused 
"by attaching other torsion tubes to form multicellular 
tubes. In the general case, tho formula for twist "becomes 



or 
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where & e is the effective shear stiffness* It should "be 
noted that longitudinals have no influence in the simple 
torsion problem. 

The derivation of these torsion formulas, which may 
■be found in a numher of standard textbooks, is based on 
the assumptions that the torques are applied as shear 
stresses distributed over the end faces according to the 
theory and that the cross sections are free to follow the 
tendency to warp that exists in most cases. In practical 
structures, it is usually not possible to comply with thos 
assumptions. The root section is usually built in more or 
less completely, and the resulting restraint on the warp- 
ing causes normal stresses, or "bending stresses, and a re- 
distribution of shearing stresses (referonco l). These 
effects disappear quite rapidly with increasing distance 
from the root and are usually negligible at a distance 
from the root equal to, or greater than, the width of the 
box. At the root, however, they may be quite appreciable. 

The commonly accepted theory of shells in combined 
bending and torsion uses the simple theories of bending 
and torsion. Corrections may therefore he necessary near 
the root to account for the effects of restraint against 
warping just mentioned, which modify the simple theory of 
torsion. Corrections may also he necessary to account for 
the effects of shear deformation, which modify the simple 
bending theory (reference 2 ) . 

An emphatic word of warning must be given relating 
to the use of the theories- of bending and torsion. These 
theories give fairly reliable results,, if they are used 
with judgment. The theory that the entire cross section 
acts as a unit naturally cannot.be expected to hold very 
well if the joints are not perfect or if - the changes in 
dimonsions and shape are too sudden. Nose covers attached 
with piano hinges and trailing edgas with their acutre an- 
gles at the tips are the most usual examples of structural 
components that cannot be expected to- be fully effective 
either in bending or in torsion. 



Sign Conventions 

External forces will be taken as positive when act- 
ing upward.- External torques will he taken as positive 
when clockwise. 



N.A.C..A. Technical Note No. 691 



7 



Bow's notation trill "be used to designate cells and 
walls. The cells will "be designated "by letters from left 
to right, starting with 11 a. u 

Shear stresses and forces in the walls of a cell will 
he taken as positive when going clockwise around the cell. 
If a wall - belongs to two cells, the sign will be estab- 
lished "by assuming the wall to "belong to the left-hand 
cell . 

Line integrations will be performed in a clockwise 
direction. It should he noted that, in a wall belonging 
to two cells, the sign of the shear must be reversed in 
the right-hand cell when performing a line integration,, 
because the arrow of the positive shear direction as es- 
tablished will oppose the sense of positive direction of 
integration. In doubtful cases, and preferahly in all 
cases, free-body diagrams should be drawn indicating the 
directions of all forces. The use of such diagrams will 
materially reduce the chances of errors^ in sign and will 
do away with the necessity of adhering rigidly to a set 
convention of signs, provided that care is taken in writ- 
ing the equations of equilibrium. -" 

Note that, in cases where the sign convention is ad- 
hered to rigidly, the basic equation (3) must be written 



* - - 2A 
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The Open Shell 

Open shells (fig. 4) can "be analyzed "by applying 
formulas (l) and (2). After the shear forces in each 
part of the cross section have "been found, the resultant 
of the internal shear forces can ho found hy ordinary 
statics. This resultant lies on a vertical line distant 
e from tho open wall of the shell. The point whore this 
resultant intersects the horizontal axis is called tho 
"shear center." The external load P must pass through 
the shear center if~thero is to he no torsion. Tho tor- 
sional stiffness and strength 6f~an open section boing 
extremely small, it is necessary to koop the external 
load vory close to tho shear center. Tho knowlodgo of 
the -shear center is therefore important for an opon sec- 
tion. 





Figure 5. 



In curved shoots with a 
constant shear intensity q., 
such as the webs of tho sec- 
tions shown in figures 4(b) 
and 4(c), it is often conven- 
ient to replace the shear 
stresses acting on the curved 
cross section by a resultant 
force. Integrating horizon- 
tal components, vertical com- 
ponents, and moments of the 
elementary shear forces, 
there is obtained (fig. 5) 
the horizontal resultant 
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H- = 0. 

the vertical resultant 

Y = gh (5) 

and the torque moment ahout any point 

T = 2qA Q (6) 

which gives as the location of the resultant force 
R = Y the distance 

. .2^ .. <?) 

h 

from the open face of the shell. In these formulas, 
A 0 is the area included "between the contour of tiye sheet 
and the open face. It should "be noted that the formulas 
do not apply to the entire sections shown in figures 4(d) 
and 4(e) subjected to "bending loads, "because the shear 
intensity would not "be constant. The formulas would ap- 
ply only to the part of the section included "between the 
two longitudinals next to the neutral axis. 

Numerical examples 1 to 5 illustrate the analysis of 
open sections. ' 
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The Two-Flange ,. Single-Cell Boz (D-Section) 

The two-flange, single-cell section (fig. 6) may "be 
considered as a combination of a beam and a torsion tube. 
The beam can take bending moments only in a plane paral- 
lel to the plane of the two flanges, so that the load P 
producing the bending must be parallel to this plane. The 
torsion tube can take care of the torsion existing if the 
load P is not applied at the shear center of the shell. 

The-total shear force acting on any cross section may 
be resolved into two components (fig. 6(b)): the shear 
force S w acting in the plane web, and the shear force 
Sjj in the nose sheet acting at the shear center of the 
nose sheet; these forces are known as to location and di- 
rection but unknown as to magnitude. There arc available 
two equations of static equilibrium to find them: 

giving S N 



and, finally 



The flange stresses are found by using equation (l ) , which 
simplifies to 



P - S ff + s N = 0 



Pd + Sjj-e = 0 



e 



2A 



a 



dh* 



P 



ht 



2A a*N 



ht 



t ff Vh 2A a ; 



(8) 
(9) 



Ay being the effective area of the flange. 





Figure 7, 
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If the angle of twist per unit length dz is desired, 
it can he found by substituting (8) and (9) into formula 
(4b) 



assuming Gg = & for both webs and denoting by p the 
perimeter, or developed length, of the nose contour. 

Prom equation (9) it can be concluded that the shear 
stress in the plane web becomes zero if P is located at 




which is the location of the shear center of the' nose web 

alone. Vice versa, if P is located at the shear center 

of the plane web, i.e., at the plane web or d = 0, the 
shear stress in the nose web becomes zero. 

Although the section shown in figure 6 is the most 
common example of the two-flange single-cell section, it 
is not the only one. Figure 7 shows another example, a 
two-spar box with three-point attachment. In this case, 
the spar attached with a single bolt cannot act as a boam, 
and the box might be termed a "rectangular D-section." 



Analysis of the D-Section by the Shear-Center Method 

If the distance d is chosen so that the angle of 
twist 0 becomes zero, the condition of bending without 
torsion is obtained, and the distance d 0 thus found lo- 
cates the shear center, or the 'elastic center, of the en- 
tire section. Letting 8=0 in equation (10), 




(10) 




2Aty 
N + P*W 



(11) 



A load P applied at d D wilT cause only bending, the 
associated shear stresses being 



N.A.G.A. Technical Note No. 691 



13 



A load P applied at a distance d can ."be replaced by an 
equal load P at d Q and a torque 

I = P(d - d 0 ) 

The shear stresses caused by the torque can be calculated 
by equation (3) and added to the stresses given by equa- 
tions (12a) and (12b). Obviously this method of analysis 
using the shear center is much more laborious than the di- 
rect method of- analysis. • ' 

Example 6 in the appendix illustrates the analysis of 
D-sections. 

The Two-Cell Torsion Tube 

Since a single-cell tube is sufficient to take tor- 
sion, a two-cell tube (fig. 8> is statically indeterminate 




Figure 8. 
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In order to analyze it, imagine the two cells to he split 
open and- unknown shears of the intensities q a = X and 

q^ = Y applied to' the cells. One equation for these two 
unknowns is furnished hy.the static equation stating that 
the sum of the torques must he zero; using equation (6), 

X X 2A a + Y X 2A b + I = 0 (13) 

An additional equation is ohtained from the condition that 
elastic continuity must he preserved-, namely, that the an- 
gle of twist of cell a relative to 'coll b must ho zero, 
or that the twist of cell a must equal that of cell b. 

e. a = e b , (14) 

Expressing— 9 a and 6^ as functions of X and Y, "by 
using oquation (4h ) , as 

6a = 7TT- / ■ •£ de and ,9-j. = A. | ds 

a 2A a *a G- D 2A-1;, ^b 5 

The second relation needed to find X and Y is ohtained 
by oquating these two expressions in conformance with (14). 

If the angle of twist "is desired, it is found by sub- 
stituting the values for X and Y in the expression for 
8 a or for 8-^ . 

Example 7 illustrates the numerical procedure. 

The Two-Flange Two-Cell Shell 

The analysis of the two-flange two=-c'ell shell in com- 
bined bending and torsion (fig. 9(a)) is closely analogous 
to the analysis of the two-cell torsion tube. Imagine the 
two cells to be split open (fig. 9(b)) and the shears of 
intensities X and Y applied. The load P located at 
d is replaced by a load P located in the plane of the 
flanges and a torque Pd. The shear intensity in the shear 
web is then 

q w = X - Y + g- (15) 

Equations (13) and (14) are again' used to find the shear 
intensities X and Y, as in tho case of the torsion 
tube; the only difference between the two cases lios in 
tho appearance of the term P/h in the web shear intonsity. 
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■It should "be noted in figure 9(a) that there is only 
a single "bolt attaching the auxiliary rear spar and that 
the flanges of the rear spar are dotted, indicating "that 
they do not enter into the calculation. 

Example 8 illustrates the analysis of a two-flange 
two-cell shell. 

■ A p 





" Analysis of the Two-Flange Two-Cell Shell 
by the Shear-Center Method 

The location .of the shear center is found as "before 
from the condition that the angle of twist must be zero. 
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Equating 6 a and €-{j each to zero, two equations are ob- 
tained instead of the single equation (14). These two 
equations together with (13) are sufficient to find the un- 
known location d 0 at which P must be placed to produce 
bending without torsion as well as the shear stresses asso- 
ciated with this special case of bending. 

After the shear center has been found and the solu- 
tions for bending only and for torque only have been com- 
pleted, any additional loading case to be investigated may 
be broken up into a combination of bending only and tor- 
sion only, as discussed for the two-flange, two-web shell. 
The analysis of any given case then consists merely in 
multiplying the stresses from the two basic solutions by 
appropriate factors and adding them. This method requires 
less numerical work than setting up and solving equations 
(13) and (14) for each case. Consequently, the shear- 
center method of analysis saves time if a sufficient num- 
ber of cases are investigated, so- that the total time 
saved on individual cases overbalances the time required 
for finding the 3hear center and making the basic solutions. 

It might be pointed cut- that the same advantages can 
"be had by using any arbitrary-load case and the pure- 
torsion case as basic cases. If the arbitrary case chosen 
as basic is for a load P x located at d x , then a load 
P 3 located at d 2 can be replaced by a load P a at dj. 
and a torque P 2 (<X a ~i. x ). The analysis of additional cases 

is therefore just as simple as if the shear-center method 
had been used. 

The shear-center method is illustrated by example 9. 

The Three-Plange Single-Cell Shell 

The single-cell shell with three flanges is of inter- 
est as the practical example of a B-section capable of 
taking bending in any plane (fig. 10). This section can 
be easily analyzed for the general case of a section with- 
out symmetry. The location of the resultant shear in each 
web is known from formula (V), and tho equilibrant of the 
load P can be resolved into three, forces alnng these 
lines by statics. If the load is parallel to the plane of 
two flanges, the third flange is unstressed, and tho shear 
intensity q is .constant for the two webs joining tho 
third- flange. The analysis is then analogous to that of a 
two-flange shell. 
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b. 




a 

Figure 10 



, The Multiflange Single-Cell Shell 

The four-flange "box (fig. 11) may "be considered from 
two points of view. If the upper cover is cut, the lower 
cover will also drop out of action. The structure is then 
the familiar two-spar wing. This structure is statically 
determinate (for vertical loads), torsion "being taken care_ 
of by one spar bending down and the other one bending upi . 
With the cover intact, the structure is statically inde- : 
terminate . 

In the commonly accepted .shell theory, the torsion 
taken by opposite bending of the spars is neglected. Tor- 
sion is assumed to be absorbed entirely by the four walls 
acting together as a torsion tube. All flanges are assumed 
to act as a unit, namely, a single beam obeying the engi- 
neering theory of bending. The shell is then, in princi- 
ple, analogous to the D-section, consisting of a combina- 
tion of a torsion tube and a beam, and is statically deter- 
minate. This conclusion remains valid if there are longi- 
tudinals attached to the cover sheets. Each longitudinal 
introduces one more unknown shear stress in the sheet and 
also one additional equation of equilibrium of forces along 
the z aii s . 

The justification for using this theory in preference 
to the one first mentioned lies in the fact that, fo.r all- 
metal stressed-skin wings, the torque taken by differen- 
tial bending of the spars is very small compared with that 
taken by the torsion tube except near the root. In the 
region of the root, corrections must be made to allow for 
this effect, as will be discussed later. 




As example, the equation for- -the simple case of figure 
11 will be developed. 

ST: P + S x - S 8 = 0 ' (a) 

ZM: -Fd + S 3 h + S a w = 0 (b) 

s s 

SL on flange li dlS + -r~ dz - — dz = 0 (c) 

n w 

S S 

SL on flange 2: dP 2 + -~ dz - dz = 0 (d) 
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Since the stresses in the two flanges are assumed to "be 
equal 

p A x dF a p A 3 

dz~ = h A~ a dz~ = h A~ (e) 

where A = k x + A a . 

Substituting (e) into (c) and (d) 

hi i+ h i -'ir- 0 (f) 

^£ + ^.% = 0 ( g ) 
hA xr h 

PA a w 

Prom (») S 3 h. = S a w - — j— 

Substituting into ("b) 

PA a w 

-pa + s s w j— + s s w = o 

p /a a s 



s = = Is ( 4 - » x) (18 > 



In tho more general case of the trapezoidal' "bos: (fig. 12(a)), 
the equations "become 

hi - h a 

ST: P + S x - S a - S 3 = 0 

EM: -Pd + S 3 h 1 1 + S a w = 0 

S, . ' S, , 

ZL X : dF, + ~ 4z t dz = 0 

x n-L w 1 

S g S a 
SL a : dP a + — ^ dz - ^— dz = 0 



In this case 

nr = ~§t~ r antt it* = tt 



dP, h,A, dP a h B A 2 _ 

1 1 P and s~~ = -4t~ P 
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so that 



and finally 



b 2 ™ - ~2j— A 2-^ 



So = P 



h s j~d + hxhgAs 
hj. + h s [w 21 



(19) 



Substituting S 2 into the preceding equation givos S 3 , 
and the first equation then yields 




figure 12. 



If a number of stringers with a total cross-sectional 
area A 3 are uniformly distributed along the width of- the 

cover sheets (fig. 12(b)), the formula for S n becomes 



S 2 = 



hi +■ h a l w 



'a 

"a kika A 2 ki (h 1 + 2h 2 )^. 3 
— + — - — + 
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J 



(19a) 



Since the analysis of a given case consists merely in 
substituting numerical values into equations (16), (17), 
and (18), no examples will be given here. If examples aro 
desired, they may bo found in reference 3, which covers in 
detail the analysis of~fche four-flange box by the shear- 
center method. - 
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The Multicellular Shell 



Inasmuch as the multiflange single- cell shell is sta 
ically determinate, the multicell shell is statically in- 
determinate. The method of 'analysis is analogous to that 
used for the two-flange, multicell shell and will "be il- 
lustrated .."by ■ the example of a four-flange » two-cell shell 
(fig. 13(a)). . . ■ ■ - 

A. 





Figure 13. 



Imagine each cell cut open and shears of intensity 
2 and Y applied. The transverse load P causes shear 
loads S x and S 3 in the -two-'spar wohs , which are pro- 
portional to the moments of inertia of these spars 
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and S 



2 




(20) 



with 



S 




I = 




2 



The equations (4"b) for angle of twist are written down, 
and the condition of continuity- 



furnishes one equation. The second equation is found, as 
"before, from the static condition that the internal torque 
must "balance the external torque. 

Taking moments about A x 



The analysis is closely analogous to the analysis of the 
two-flange, two-cell soction previously discussed. The 
only difference is that the shear intensity in the rear 



Example 10 illustrates the analysis of a four-flange, 
two-cell section "by the direct method. Example 11 illus- 
trates the method of finding the shear center and the shear 
stresses associated with bonding. 

The procedure for more complicated casos (fig. 14) is 
merely a simple extension of th©- procedure discussed, so 
that no example will "be required. 



9 a = e b 



-Pd + S s w + 2A a X + 2A- 0 Y = 0 



(21) 



spar is now Y + r— instead of Y. 

"a 



A 




figure 14. 
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Corrections to Simple Theories of 
Bending and Torsion 

As mentioned before, the simple theories of "bending 
and of torsion used thus far may require corrections. 
These corrections are important only in the inboard region 
near the root, if they are important at all. Whenever they 
are to be made, it is very advisable to separate bending 
from torsion at the outset. Such a procedure will make 
the calculations much clearer and will materially reduce 
the danger of committing errors in sign. 

The flange material in box beams is usually distrib- 
uted across the cover in the form of individual stringers 
or corrugated sheet. The bending action of such beams 
differs from that assumed by the simple theory of bending, 
because the sheet deforms under the shear stresses imposed 
on it. The analysis of bending action under such circum- 
stances is discussed in reference 2 and no discussion will 
be given here. 

In tubes subjected to torsion, the cross sections 
usually have a tendency to warp out of their original 
planes. If this warping is prevented by attaching the 
tube to a rigid support, or by conditions of symmetry in 
the middle of the span, then longitudinal (normal) stresses 
will arise, and the shear stresses will be redistributed. 
Reference 4 gives a method of calculating the effects for 
cross sections of arbitrary shape, but the method is of 
limited usefulness. Methods of analyzing rectangular tubes 
have been developed by a number of writers; th.e most im- 
portant reports dealing with the methods are summarized in 
reference 1, This reference also gives what appears to be 
the only published experimental data. They do not agree 
very well with the theory; fortunately, the effects of end 
restraint are small in most practical cases, so that they 
need not be very accurately calculated. 
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For a rectangular ' tube symmetrical about both axes, 
such as shown in figure 15, the normal forces on the cor- 
ner flanges caused by complete restraint may be calculated 
by the formula • 




(22) 



(reference 1, equation (3c) with G-/E = 0.4). The sign of 
the stresses is determined from- the rule that the walls 
with the smaller- ratio of width to thickness (in general, 
the vertical walls) act like independent spars, absorbing 
tho torque by bending in opposite directions. 




<f b 3» 



Figure 15. 



Tho effect of tho end restraint on the shoar intensi- 
ties is written most conreniontly in the form 
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Ag. = ± 



T 

2A 



Jb_ 



t„ 



(23) 



obtained "by using oquations (3c), (9), and (10) of refer- 
ence 1; Aq. is the correction to "be applied to the shear 
intensity calculated on the assumption of no restraint, 
namely, 




The negative . sign in (23) is used for the walls with the 
larger ratio of width to thickness, in general the hori- 
zontal walls. 

In actual cases, the "box will seldom "be symmetrical 
ahout "both axes, as assumed in the derivation of formulas 
(22) and (23). The simplest procedure in such a case will 
"be to use average values for "b , c, t^ , and t c . This 

procedure is somewhat unconservative , "but formulas (22) 
and (2*) are "basically conservative "because they assume 
infinitely closely spaced "bulkheads. Furthermore, except 
in such cases as wings continuous across the center line 
of the airplane, the root section will not • "be rigidly 
"built in, "because there will "bo play in the fittings and 
olastic yielding in the fittings and in the contor sec- 
tion. 

The simple procedure outlined hero and used in ox- 
ample 12 may, of courso, he insufficient in some casosj 
a moro detailed treatment, however, is "beyond tho scope 
of this paper. 



Langloy Memorial Aeronautical Lahoratory, 

National Advisory Oommittoo for Aeronautics, 
Langley Field, Ya. , January 20, 1939. 
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APPENDIX 
Exampl-e 1 



Pind the shear intensity and the shear center of the 
section shown in figure 16 for a vertical load P. 

I = 2Az S = 2xi-x4 = 2in. 4 

ft = Ax = | >! 2 = | in. 3 

By formula (2), the shear intensity is 

q = --g-A = J P Id. /in. 

In this simple case, the correct direction of the arrow can 
he found by inspection, and the sign convention is not needod. 



The force H in the leg. (fig.' 16(b)) is 
H - qw - J P X 3 = J P 



H X 4 = 3: X P 



x = 3 in. 



or the location e of the shear center e=x+w=6 in. = 2w, 
which agrees with formula (7). 



h=4 



t=.065 ," 1 
o 



. w =3 J 

(a) 



A 



Pigure 16. 



H 



3C S : 



(b) 
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Example 2 

Find the shear intensity and the shear center for the 
section shown in figure 17. 

Bow's notation is used as indicated: 

I = 2 x J x 4 + 2 X J X 4 = 4 in. 4 

^ac - \ * 2 = I 

^ah - J x 2 + i X 2 55 1 
The shear intensities are therefore 

n _ P x j _ 1 p n P X 1 

4 ac ~ 4 " 8 Vb ~ 4 

The horizontal force H is 



H = a x i = 5 P X 3 4 P 

ac .8 a r 

Taking moments about the lower left corner of the channel 
as "before , - - . 



Px = Hh 



s. P x 4 
8 



x = | in. 



e = I + 3 = 4 i ii. = |w 



t=.065 



•7V 



h=4 



w=3 

(a) 



Figure 17. 



Lb) 
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Example 3 

Pind the shear intensities and the shear center of 
the section shown in figure 18. 

I = 6 in. 4 

Qad = | ln - 3 <laa = * x ife 1^. /in. 

^ac = 1 <*ac = P * | "-/in. 

Q a -b = li in. 3 qL ab P x | lb. /in. 

H = q , X 1.5 + d. X 1.5 = t% X 1.5 + | X 1.5 = | P 

"ad ^ae 12 b 8 

Px=Eh=£px4=|P 
8 2 

x = § in. 

as in the preceding case. The relation e = ^ tr holds 

for all channels, if the effective material is uniformly- 
distributed along the legs of the channel. 
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Example 4 

Find the shear intensity and the shear center for the 
section shown in figure 19. 

a 4 
1 = 2 A p E in. 

Q = XjR in. 3 

q. =• P — — -s = -"- lh./in. 
2A p E 2R ' 

Taking moments a "bout the center 
of the circle 



Pe = i x uE x a = — x tiI x 1 



- ttPR 
2 



e = § E in. 




Figure li). 



Example 5 

Find the shear center of the section shown in figure 
20. In this case, all the sheet is effective in bending. 

T 3 4 

I = 4 ttE t in. 

2 i 

Q, e = / EtdGxR cos 6=R t sin 9 in. 3 
o 

_ R a t sin 9 _ 2 sin 6 ,, /. 
q.o = Px- 5 =P — — - — lo/ in. 



£ irE 3 t 



Taking moments about the center of' 
the circle ■ 

Pe = f 17 P — ai g 6 ' Ed 6 x E 




= — R in. 



Figure 20. 
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Example -6 



Given the D-soction 
shown in figure 21. 

(a) By diroct analysis, 
find tho stresses in the 
section and tho anglo of 
twist, assuming that no 
"buckling occurs. 

("b) Find the shear cen- 
ter of the section, and 
make the analysis "by the 
shear-center method. 

(c) Find the changes 
caused "by the flat sheet 
developing a full diag- 
onal-tension field. 

( a ) direct anal v- 
s.L£-- Find first the lo- 
cation o of tho shear 
center of the nose by 
formula (?) 




t w =.032 



P=5,000 



tjj=.064 



2A, 



2 5 b. 1 



= - B = 15.71 in. 

2 . 



U N ~ h ~ 2B 
Taking moments about the plane web 

S N X 15.71 = Pd = 5,000 X 5 = 25,000 in. -lb, 



S N = 1,591 lb. 



which gives 



Stf = 6,591 lb. 
The shear stresses are therefore 

= 1,244 lb./ sq.in. 



_ 1591 

N ~ 0.064 X 20 



6591 = 10,300 lb./sq.in. 



0.032..X 2 0 

The angle of- twist is obtained from the "basic formula (4b ) 
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9 * 2^ $ Tds 



£l,244 X tt X 10 + 10,300 X 20 J 



2 X g X 100 G- 
= | X 780 

With ft = 4 x 10 s , this value becomes 

3 = 195 x 10~ 6 ■ radian per inch length, 

(b) S hear-c e nter analysis .- In this simple case, the 
location of the shear center can he found from formula (11) 

2 X 3 X 100 X 0.064 

d_ = - — ; » - 8.80 in. 

° 20 X 0.064 + tt X 10 X 0.032 

In order to illustrate the procedure used in general cases, 
the solution will be carried through starting from funda- 
mental principles. 

The load P is assumed to act at an unknown distance 
d. There are then three unknowns: Tjj , T^, and &. To 

find these unknowns, there are available the static equation 
EV = 0 and EM = 0 and the elastic equation 9=0. 

ST = q H h - a, w h + P = 0 

2M (about the web) = q. N he - Pd Q = 0 

6 = j~ [j^ x P + x hj = 0 (by equation (4a)) 
Numerically 

q.jj X 20 - q. w X 20 + 5,000 = 0 

q. K X 2 0. X 15.71 - 5,000 d D = 0 

96 - i x hrri [rMi 10 + oJfe x g °] - 0 

These three equations are solved and yield 

q. N = - 140.1 lb./in.; q. ff = 110.0 lb. /in. 

d D = - 8.80 in (a) 
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The shear intensities just obtained are' those associated, 
with "bending caused by a load P applied at the shear cen- 
ter. The actual load is applied at d 5, so that there 
is a to roue 

T = — 5,000 (5 + 8.80) = - 69,000 in. -lb. 

giving a shear intensity 

» - * f gg?.i - F 9 - 8 ' c *> 

The total shear- intensities are then • 
q.jj = - 140.1 + 219.8 = 79.7 

and T N = q7064 = 1,245 113 '/so., in. 

ol w = 110,0 + 219.8 = 329.8 
and 1" w = q = 10/300 lb./sq.in; 

(c) Chang e s caused by "buckling in flat we b.~ Tfhen the 
flat web is allowed to.buqklo into a full diagonal-tension , 
field, the effective thickness becomes 

t a = | t = | X 0.032 = 0.020 in. 
° 8 8 

In order to evaluate the angle of twist 8, it is necessary 
to obtain an effective shear stress 

T e = f T = f x 10,300 = 16,500 Ib./sq.in. 

Substituting this value into the expression for 9 

6 = Yoo^'g j^ 1 * 244 it X 10 + 16,500 X 20 J 

- — x 1,174 = 293.5 x 10~" 6 radian per inch 

The changed location d D of the elastic center is obtained 
by substituting t e instead of t into formula (11 ) 

2 x J x 100 x 0.064 - 

d A = : = - 10.55 in. 

, 20 X 0.064 + tt X 10-X 0.02p . 
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Example 7 ' - 

Find the shear ■ stresses and the angle of twist in the 
torsion tube shown in figure 22. Assume G 0 = G- for all 
walls. 



t=.020 t=.0V3 
\ 




n 1= 24 /V=392 



i-s— 



23 



c __t=.051 
^ A b =990 

t=.030 



t=.036 



h P =20 

A 



44 



(a) 




Figure 22. 



It will he necessary to sot up expressions for tho an- 
gles of twist 6 in terms of the shear intensities X and 
Y by using formula (4a). As a preliminary stop, tho auxil- 
iary parameters 

. a = / da 

■> t - 

will ho computed, so that formula -(4a) will ho used in the 
form 

8G = a* 



34 K.A.C.A. Technical Note No. 691 



Bow's notation is used as indicated. 



44 - 602.5 



4 Dd ~ 0.073 



a be ~ o!o36 " 555 



01 0.030 



ha 0.051 



(Note that a-^ = a a - fa . ) 

The expressions for the angles of twist are 



6„G = 



fx x a a . n + (X - Y) a a J 



2>TA a {_ ac a *_ 

= 3.85 X - 0.600 Y .. 



9 0 - — 
a 2A 



= 2^~q j^Y X 602.5 + Y x 555+ Y X 1,467 + (Y-X) 470 J 

= 1.562 Y - 0.2376 X 
Equating 8 & = 9^, obtain 

4.088 X = 2^162 Y (a) 
The equation of moment equilibrium (13) is 

X x 2 x 392 + Y X 2 X 990 - 250,000 = 0 (b) 
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The solution of these two equations is 

X = 55.3 lb. /in. Y = 104.5 Id. /in. 

The shear stresses are therefore 

T = S 55,3 = 2,765 lb./sq. in. 
as t 0.020 * 

T oe = f = £M = 2 ' 900 



T _ I _ iPAii = 3,480 lb./sq. in. 
01 t 0.030 

X_ = _Y = _ 49.2 = _ 965 i^./ S(1 .in. 



"f — _ ii ^ ■ * - 

ab t 0.051 

The angle of twist is obtained by substituting into the 
expression for 6 & 

ft = 6 a = i |"3.85 X - 0.600 y] = X 150.3 

With G- = 4 x 10 6 this value becomes 

G = 37.6 x 10~ 6 radian per inch length, 
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Example 8 

Find- the shear stresses in. the section shown in fig- 
ure 23. This section is identical with the one used in 
example 7, except that two flanges have "been added to take 
care of beam action. 

^ P=5,000 




(a) 




The parameters a can he taken from the preceding 
example . 

The shear intensity caused by the load P acting on 
the shear web is 

The expressions for 9 are now written exactly as in ex- 
ample 7 except for the addition of . 
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8 a, & = 2~x"392 [ X * 2 ' 550 + (X - T + 208.2) X 4?oJ 
= 3.85 X - 0.600 Y + 124.96 . 

6-tjG' = 2~^q |\x602.5 + Yx555 + YX1467 + (Y-X-208.2) X 47oJ 

= 1.562 Y - 0.2376 X - 49.5 
Equating 6 & & to. G^G- gives 

4.088 X - 2.162 Y + 174.46 = 0 (a) 

The equation of moment equilibrium is taken around the 
shear web, to eliminate one term ■ 

X x 2 X 392 + Y X 2 X.990 - 5,000 X 46 = 0 ' (b) 

Solving equations (a) and (b), obtain 

X m 15.5 Y = 110.1 

The shear stresses are therefore 

T o~ = TT 1 ^ = 775 Ib./sq.in. 
ac 0.020 . 

= 0T073 = lj51 ° • 

T * e = oTosi 3 3 ' 060 ll5 -/sq.in. . -■- 

= So30 = 3,670 Ib./sq.in. 

■ T at = .15^-5_- i L S:, + 208,2 = 2,226 lb./sq.in. 



Example 9 

For the section used in example 8, find the shear 
center, and analyse the load 'case of example 8 by the shear- 
center method. 

Tiro equations are obtained by equating to zero the ex- 
pressions for 8 0 and 8>.," which are taken from example 
8. 
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" 3.85 X - 0.600 Y +124.96 = 0 

-0.2376 X + 1.562 Y - 49.5 =0 

The solution of these equations gives the shear intensi- _ 
ties associated with torsion-free bending. 

X = -28.2 lb. /in. Y = 27.42 lb. /in. 

The distance d Q of the shear center from the shear web 
is obtained by writing EM about the shear web 

-5,000 d - 28.2 X 2 X 392 + 27.42 X 2- X 990 = 0 

d Q = 6.44 in. 

A load P located at d = 46 inches- will therefore cause 
a torque 

-P(d - d Q ) = -5,000(46 - 6.44) = -197,800 in. -lb. 

The stresses Tg due to bonding are obtained from X and 
Y as before. 

The stresses due to the torque of_-197,800 in. -lb. 

are obtained by multiplying the stresses from example 7 by 

197800, _ n 7Q1 
250000 " °- 791 * 

The final stresses T are obtained by superposition as 
shown in tho following table. 



Wall 


T B 

(lb./ sq. in, ) 


T T 

(lb./sq. in. ) 


T 

(lb . / sq. in. ) 


ac 


-1,410 


2,190 


780 


bd 


376 


1,132 


1 ,508 


be 


762 


2,297 


3,059 


bf 


914 


2,757 


3 ,671 


ab 


2 ,990 


-764 


2 ,226 
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Example 10 

For the section shown in figure 24, find, the shear 
stresses. The section is identical with that shown in 
figure 22 except for the flanges. The load is assumed to 
he perpendicular to the neutral "axis. The inclination of 
the shear we"bs is neglected. 



a P=5,000 

' A 2 =1.65 
*1 





d 




Figure 24. 



Imagine the cover to he slotted in "both cells as in- 
dicated," leaving a structure consisting of two spars. The 
vortical shear is divided hetween these two spars in the 
ratio of their moments of inertia. 

I x = | X 1.85 X 24 s = 533 in. 4 
I 3 = i 'x "1.65 X 20 3 "= 330 . in. 4 



I 



= I 1 + I s = 863 in. 4 
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There f oro 



, - P X = 5,000 X Iff = 3,090 lb. 
■ 1 obo 

S x 3090 



I, »jj » -34- = 128.8 lb. /in. . 
S a = P X ^ = 5,000 X HI = 1,910 lb. 

Prococding as in the previous cases, write the expressions 
for 9 a and 8^. 



9 a - b 



1 2,550 X + 470 (X - t'+ 128.8) 

X <52d 1 J 



= 3.85 X - 0.600 Y + 77.3 

6-b - 2x990 [ 602 - 5 *+ 555(1+95.5) + 1,467 Y + 470 (Y-X-128 . 8 ) J 

= 1.562 Y - 0.2376 X ~ 3.8 
Equating 6 & to 6^, obtain 

4.09 X - 2.162 Y +81.1 = 0 
The equation of moment equilibrium gives 

-5,000 X 46 + 1,910 X 44 + 2 X 392 X + 2 X 990 X Y = 0 
Solving these two • equat ion s , obtain 

X = 15.9 lb. /in. Y = 67.5 lb. /in. 

The shear stresses are therefore 



X 



ac 0.020 



= 795 lb./sq.in. 



t ^a = n = 925 lb./sq.in. 

ba 0.073 ' - 
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- ^T^r = 4,530 "•/■i-. in - 

T ^ = od)30 = 2 ,250 lh./sq.in. 
X.+ g. x - Y . , ■ 



Example 11 

Find the shear center, and the shear stresses associ- 
ated with torsion-free tending, for the section shown in 
figure 24. 

Take the expressions for 6^ and _ 6-5 from the pre- 
ceding example and equate each one 'to zero. 

3 ...85 X - 0.600 Y + 77.3 ^.0 

• * •■ 

-0.2376 X + 1.562 Y - 3.8 .= 0 

Spiring 

X '= -20.2"lh./in. Y = -0.635 lh./in. 

The shear stresses are therefore 

~- ■ -T ac 57§f§ = -1.010 Ih./sq.in.' 

T hd £^4Sf = ~ 9 lh./sq..in» 

Da ' 0.073 

-0.63 5 + 95 .5 „ 

T be = d"7o36 = 2,640 113 -/Esq.- in - 

T *f = -§7o§0 = " 21 ^-/Bd.in. . 

- ' ' T a . = -20^ 2 + 128,8 + 0.6 55 = 2 14Q n,./ sa _. in . 
. a ° 0.051 ' 

Lea-ring the location-' d of. the load P undetermined, 
write the moment equation 

EM = -5,000d + 1,910X44 - 2x392x20.2 - 2x990X0.635 = 0 
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which gives as location of the shear center 

& = 13.39 in. "behind front shear web 

Example 12 

For the section analyzed in examples 10 and 11, find 
the stresses if the section is a root section that is rig- 
idly built in. The length L of the beam is 200 inches,* 
the load P is applied at the tip. 

The first step is to separate the load on the entire 
section into bending moment and torque. The "bending mo- 
ment is 

M = PL = 5,000 X 200 = 1,000,000 in. -lb. 

According to example 9, the shear center is located at 
d 0 = 6.44 inches, and the torque is 

T = P(d -'d 0 ) = -197,800 in. -lb. 

The effects" of restraint against warping will be calcu- 
lated under the assumption that only the approximately 
rectangular cell "b "between the four main fittings is re- 
strained against warping and that the nose part has no 
influence on these warping stresses. 

According to example 7, a torque of 250,000 in. -lb. 
creates a shoar intensity Y = 104.5 lb. /in. in coll "b . 
The existing torquo ofi'197,800. in. -lb. will thorefore giro 
a shoar intonsity 

q,, = 104.5 X i|J|§§- = 82.8 'lb. /in. 
*D 250000 < 

The torque carried by cell "b .is therefore (approximately) 

Tfc = 82.8 X 2 x 990 = 164,000 in.-lh. 
With the average- values 

•b = 44 in. t^ = 0.075 + 0.030 = 0<0515 in# 

o = 24_|_20 m 22 in> ^ m 0.051 | 0.036 = 0t0435 in< 
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. Aj. = If 85 + 1.65 „ 1.75 . S q..in.. 

The normal force on the flange due to torque "becomes, by- 
formula (22) 



X = 0.56 164000_ (855-506) / ■ 1,75 e ~ = 1,190 lb. 
44 X 22 . V 855 + 506 

and the correction for shear intensity "becomes, "by formula 
(23) 

An = 164000 (855 - 606V 21 . 75 'lb./in. 

H 2 X 44 X 22 (855 + 506) 

The "bending stresses due to the "bending moment are, in the 
front flanges, 

a-- = ± 1000000_x_12 = ±1 i5 >9 oo lb./sq.in. 
±JS 863 

and, in the rear flanges, 

- = ± 100QQ00 X 1 0 ^ ±11,580' lb. /so., in. 

2B • 863 ■ 1 

the upper sign applying to the upper flange in each case. 
The bending stresses duo to torque are 



and 



" k " lift ^^S lb./sq.in. 

a 2i -17 ■ - I7§§ =r±721 "-/^a.in. 



the upper sign applying to tho upper flanges. The final 
stresses aro therefore. • " ." : " "~ . 

■cr 1IJ = -13,900 + 643 = -13 ,25.7' lb-./ sq. in* 
cr 1L = 13,900 - 643 = 13,257 lb./sq.in. 
a 2U = -11,580 - 721 = -12,301- lb./sq.in. 
CT 2L ~ 11 » 580 + 721 = 12 * 301 3 -' b »/sq.-in-- 
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The shear stresses for the section free to warp are ob- 
tained from example 10. Superposing the corrections 
Aq. 

At = gives the final shear stresses 

T ac = 795 + 0 ~ 795 lh./sq..in. 

T od = 925 - §7$H = 627 "-/Bq.ln. 
T^ e = 4,530 +.§"§§ = 5,134 lh./sa.in. 
T M = 2,250 - = 1,525 lb./sq.in. 

T ba = -1,513 + fp^rf- = -1,087 lb./sq.in. 

The corrections influence the design of the rear spar more 
critically than 'the front spar. The correction on the 
flange stress is somewhat over 6 percent; on the web shear 
stress, it is somewhat over 13 percent.. An error of- 25 
percent on tho correction would therefore cause an error 
of 1-1/2 percent on the flange stress and an error of 3 
percent on tho web shear stress. 
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